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ABSTRACT 





Seismic analysis of earth and rockfill dams is generally done in two ways: quasi-static and dynamic. However, 
a quasi-static method with easy application and simple assumptions may lead to unsafe and uneconomical 
results. In the present study, two static and dynamic analyzes have been used nonlinearly using the Rayleigh 
Damping rule to calculate the stress and strain of Azadi Dam in the stages of the end of construction and 
steady-state seepage. Also, in numerical analysis, Abaqus software and a simple elastoplastic behavior 
model based on the Mohr-Coulomb criterion have been used. The results show that in both quasi-static and 
dynamic seismic analysis, the highest strain of the Azadi Dam core occurred at the upper levels of the core 
and the highest stress occurred at the level of the core floor. The stress in the dynamic state is higher than 
the quasi-static one in the directions 0xx 49%, Oxy 30%, and Oyy 28%, respectively. Also, the maximum shell 
stress at 1255 m, 1275 m, and 1300 m levels is 29%, 68%, and 72% higher than the core, 
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1. INTRODUCTION 


eismic design and analysis of earth and rockfill 

dams are done by two methods, quasi-static and 

dynamic. The method of dynamic analysis is 

mainly based on stress analysis and displacement, 
which is usually done with the help of finite element 
methods. This method is commonly used to analyze the 
stability of large dams in the study phase. Lack of accurate 
software for dynamic analysis of earth dams, the limited 
number of experts aware of dynamic analysis, the 
complexity of dynamic analysis method, expensive tests 
for determining dynamic soil properties, frequency, and 
ease of analysis with quasi-static software are the reasons 
for widespread use of the quasi-static method. Due to these 
cases, determining the accuracy of the quasi-static method 
and creating a relationship between the solutions of the two 
quasi-static and dynamic methods 1s of interest to earth and 
gravel dam design engineers. Today, the development of 
finite element and finite difference software has made it 
possible to use dynamic analysis as well as quasi-static 
analysis. Ambraseys and Sarma, 1967examined the 


response of earth dams to several earthquakes. They 
calculated the time history and distribution of earthquake 
acceleration in the dam body [1]. (Sarma, 1975) developed 
diagrams for calculating the critical horizontal acceleration 
in which the critical horizontal acceleration is_ the 
acceleration that can bring the soil mass limited to a 
landslide level into equilibrium [2]. (Wang et al., 2006) 
introduced a new model in FLAC software and 
dynamically analyzed several earth dams in the effective 
stress space. They compared the actual deformations of the 
dams with the estimated values with different models [3]. 
(Tsai et al., 2006) by studying the dynamic response of the 
Pao-Shan dam, studied the effect of core dimensions on 
the potential of earth dam response as well as the effect of 
core width and height ratio and dam length and height ratio 
at the first natural frequency [4]. (Tsompanakis et al., 
2009) Using a neural network, evaluated the dynamic 
response of the sample embankment (laboratory) using the 
finite element method. Considering the nonlinear behavior 
of soil materials, he concluded that the magnification 
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module would shrink as the maximum earthquake 
acceleration increases and the materials enter the nonlinear 
section [5]. (Elia et al., 2011), investigated the seismic and 
aftershock behavior of the Marana homogeneous dam in 
Italy [6]. (Mukherjee, 2013) The basic concepts of 
different methods of seismic stability analysis of earth 
dams and salient features, advantages, and limitations of 
each. To realistically predict the earth dam response during 
an earthquake, the factors of nonlinear soil elastic 
behavior, the dependence of the enclosed soil pressure on 
its stiffness, the geometry of the valley, and the 
intersection of the dam with the alluvium must be carefully 
considered [7]. (Huang, 2014) analyzed the seismic 
response of earth dams with stabilizing materials 
(materials with low strength control, CLSM) using the 
finite element method. Their results show that the use of 
CLSM is suitable for stabilizing embankments against 
seismic excitation [8]. According to studies (Panulinova 
and Harabinova, 2014), the stability of earth dams against 
landslides or seismic effects should be designed so that the 


2. MATERIALS AND METHODS 


embankment is not destroyed due to changes in soil 
properties or external influences and remains stable [9]. 
(Bandini et al., 2015) presented a limit equilibrium model 
in which changes in block geometry and changes in shear 
strength due to slip are considered. He compared the 
results of observational models with the results of 
numerical analysis [10]. In all these comparisons, the 
observed behavior is consistent with the predicted one, 
which indicates the need to consider the block geometry 
change and the shear strength due to shear in the 
calculations. In this study, quasi-static and dynamic 
analysis of stress-strain in the Azadi earth dam after the 
end of the construction phase and in the steady-state 
seepage using Abaqus software and nonlinear analysis 
have been investigated and compared. Finally, the degree 
of conformity with the percentages presented in the stress- 
strain analysis at different levels of the body and core of 
Azadi Dam, as well as the difference between quasi-static 
and dynamic methods, are shown. 


2.1. GEOGRAPHICAL LOCATION AND GEOLOGY OF AZADI DAM 


Azadi Reservoir Dam is located in Iran and Kermanshah 
province, 500 m downstream of Shahgozar Bridge, and 
about 90 km from Javanrood city in the coordinates of 
46°21 east longitude, and 34°33 north latitude on Zamkan 
River. The catchment area of this river up to the axis of 
Azadi Dam is 1054 km2. Access to Azadi Reservoir Dam 


is possible through Kermanshah Road, Kuzran-Shahgozar 
Bridge. Azadi Dam is a rockfill type with a clay core with 
a height of 64 m from the foundation to 1312 (masl) meters 
above sea level and a crest length of 737 m. The volume 
of the dam reservoir at the normal level is 57.47 million 
m3 and the useful volume is 50 million m3 (Figure 1). 





Figure 1. Azadi earth dam, Kermanshah, Iran 


Azadi Reservoir Dam site consists of two rocky sections 
of Chilean Amiran sediments and marl limestone of Gurpi 
Formation and the alluvial-debris sediment unit of the 
present period in the right-axis ridge and under the 
overflow. In terms of geostructure, the area of Azadi Dam 
and related facilities belong to the folded Zagros structural 
unit in East Lorestan. This structural unit, like other folded 
Zagros regions in the south and southwest of Iran, has a 
stressful past in terms of tectonic activity. Of course, in this 
area, folding and faulting of formations are normal and 
natural (Abdan Faraz Consulting Engineers). Azadi 


Reservoir Dam site is composed of two rock sections: 
Amiran Formation shale sediments - Gurpi Formation 
marl limestone and alluvial-debris sediments of the present 
era in the right-axis and below the overflow. In terms of 
geostructure, the area of Azadi Dam and related facilities 
belong to the Zagros fold and thrust belt structural unit 
(Zagros FTB) in East Lorestan. This structural unit, like 
other Zagros FTB areas in the south and southwest of Iran, 
has had many tectonic activities in the past. Of course, in 
this area, the folds and faults of the formations are normal 
(Abdan Faraz Consulting Engineers) [11]. 


2.2. GOVERNING EQUATIONS FOR STRUCTURAL DYNAMICS 


By discretizing the dynamic equation of the structure and 
considering the applied forces of the earthquake in the time 
domain, and using the finite element approach, the 


[M]{U%C..) J+[C]{U%. + [Kk] {UJ={F_I }-[M]{U_g }+[Q]{P} 


dynamic equation governing the dam and foundation will 
be written in matrix form (1): 


(1) 
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[M], [C], and [K] are the matrices of mass, damping, and 
stiffness of the structure, respectively. {U},{U%.},{U%C.) 
}, {F_1 } and {U_g } are the relocation, velocity, structural 


acceleration, body forces, and earthquake acceleration, 
respectively. 


2.2.1. Complete Elastoplastic Analysis of Embankment Assuming Mohr-Coulomb Theory 


In an elastic-plastic analysis (complete plastic), the 
beginning of the stress-strain curve 1s linear, and its plastic 
range is linear. A yield function must be defined to 
evaluate whether the point has reached the plastic limit. 
Yield Criterion is usually expressed in terms of principal 
F=f(o_1,o_2,0_ 3,n_1,n 2,n 3) 


ni show the direction of the main stresses oi. If the 
materials 


F=f(o_1,o_ 2,0 3) 


If the stress field is such that F (o) <0 is the behavior of the 
elastic material, and as soon as the yield point 0 = F (0) is 
reached, the plastic behavior of the material begins. In the 


de= [(de]] *e+ [del] Ap 


For plastic strains, the law of flow is determined. The law 
of flow assumes that the plastic strain is perpendicular 


(de) “p=A (0 f(0))/( 0 0") 
Where A is scalar and f (o) is a level of stress function. If f 


(o) is the same as the yield function, the related flow law 
holds. Otherwise, the law of flow will be unrelated, in 


(de) “p= 0 g/( 0 0”) 


A: is called the plastic coefficient, which in the elastic 
condition has a value of zero and in the plastic condition 


o*'=[D“%e-o/d D’e 0 g/(d 0%') (0 fT) 0 0 ) D*eJE*0 


d=(0 fT)\(0 0") Dre Ag/(d 0”) 


If the soil behavior is elastic, a is zero, and otherwise, a 
equals one. Also, f is the yield function and g is the plastic 
potential level. If the Mohr-Coulomb Yield Criterion 1s, 


f 1=12 lo 2+ bo 3”'4 I+ 1/2 (o_2”'+o_3”' )sing-c.cosp=0 


f 2=12 lo 3+ bo I” I+1/2 (o_3”%'+o_ 1” )sing-c.cos@=0 


f 3=12 lo 1" bo 2” I+1/2 (o_1“'+o_2”' )sing-c.cosp=0 


The main parameters representing the Yield Criterion are 
the internal friction angle (@) and soil cohesion (c), 
respectively. The shape of the function is conical in that 
the points inside it show the elastic range and the border 
g_1=1.2lo_2" 4-- 0_3”'4 14+1.2(0_2”'+0_3” )sin 
g_2=1.2l0_3” 4--o_1"414+1.2(0_3”'+0_1” )sinW 


g 3=1.2lo_1" 4-- 0_2"'414+1.2(0_1"'+0_2” )sinW 


stresses or stress tensor variables. The onset of the 
condition is determined by the surrender criteria. The 
general form of the Yield Criterion can be given as 
Equation (2). 


(2) 


are the same, the Yield Criterion becomes a simple 
equation (3): 


(3) 


complete elastic model, the strain diagram consists of two 
components, elastic and plastic (Equation 4): 


(4) 


to a plane. This law is defined as Equation (5): 


(5) 


which case, in addition to defining the yield function, a 
new function [g (o)] will be defined, on which the plastic 
strain diagram will be perpendicular (Equation 6): 


(6) 


will have a value greater than zero. The general 
relationship between the effective stress diagram and the 
strain diagram can be expressed as Equations (7) and (8): 


(7) 


(3) 


the Yield Criterion is defined as relations (9), (10), and 
(11): 


(9) 
(10) 


(11) 


points show the plastic threshold. Since there is no related 
flow law in the Mohr-Coulomb Yield Criterion, the g 
function for the model is defined as a relation (12), (13) 
and (14): 

(12) 

(13) 


(14) 


le 
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Parameter ‘¥ is used to model the volumetric strains of 
plastic in soils that increase in volume during cutting. Also, 
in the presence of cohesion, the Mohr-Coulomb model 
allows the element to be stretched, but in the modified 


f 4=0_1%'-o _t0 
f 5=0_2”'-o t20 
f 6=0_3%'-o t20 


At new levels, it is assumed that the related law is in place. 
If the stress range is within the yield function, the body 
behavior will be a function of the hook linear model. 
According to what has been said, in this model, the stress- 
strain relationship is defined by defining 5 parameters that 
can be achieved by known and common experiments in 


Mohr-Coulomb Yield Model used in the program, the 
points under tension can be eliminated by defining 


complementary functions. These functions are defined as 
relations (15), (16), and (17): 


(15) 
(16) 
(17) 


soil. These parameters are soil shear modulus, Poisson's 
ratio, friction angle, cohesion, and expansion angle, which 
are formulated in equilibrium and compatibility equations 
in each of the elements by assuming planar strains and are 
determined by gradually applying loads and comparing 
them with yielding levels. (Zienkiewicz, et al., 1977) [11]. 


2.3. MODELING AZADI DAM IN ABACUS SOFTWARE 


Abaqus is a set of highly powerful finite element modeling 
programs capable of solving simple to complex linear 
analysis and nonlinear modeling problems. In _ the 
nonlinear analysis, Abaqus automatically selects the 
values of the convergence tolerances and also adjusts their 
values during the analysis to obtain the correct answer. As 
a result, the user rarely has to specify the values of the 
numerical solution control parameters. It also supports 
Python open-source programming language for 





programming within the software. The ability to write 
scripts in software doubles its modeling capabilities. In 
this research, Abaqus has been used to calculate the stress 
and strain pressure, assuming the flat strain behavior in 
Azadi Dam. For this purpose, the largest section of the dam 
has been modeled using Abaqus software and analyzed 
with eight-node elements [12]. Figure 2 shows the 
modeling and meshing of the Azadi Dam. 
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Figure 2. Modeling and meshing of Azadi Dam in Abacus software 


For dynamic analysis, it is first necessary to perform quasi- 
static analysis and after equilibrium, dynamic analysis is 
started. The damping used in the dynamic analysis of the 
Azadi Dam is of the Rayleigh Damping type. Rayleigh 
Damping is the most common type of mechanical damping 
used in dynamic analysis. Rayleigh Damping 1s generally 
used in time-dependent applications to provide attenuation 
that is almost independent of frequency. Damping 
percentage is considered equal to 1% due to the 
elastoplasticity of the behavioral model of the materials. In 
behavioral models that allow the soil to enter the plastic 
part (Mohr-Coulomb), considering the energy dissipation 
capability in the model, it is reasonable to include damping 
between 0 and 1%. In fact, for most dynamic analyzes that 
involve large strain conditions, only a small percentage of 
damping is required. To evaluate the performance and 
seismic design of dams against earthquakes, the force 
caused by the earthquake should be suitably applied to the 
dam structure and the seismic responses of the dam should 
be calculated by performing nonlinear analysis. Since the 
Azadi Dam site is located on the Shale rock foundation, 
the location of the accelerometer must be consistent with 
the geological conditions of the site. Therefore, for the 
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dynamic analysis of Azadi Dam, the accelerometers of 
earthquakes have been selected that have been recorded on 
rocks or rocks with a shear velocity of less than 760 m/sec. 
It should be noted that the accelerometers have been 
selected based on the type of soil at the station (soil II). For 
this purpose, the soil of the stations in question has been 
determined based on geophysical methods. Based on 
seismicity studies in the area of the Azadi Dam 
construction site, the values of design basis seismicity 
parameters (DBL), design top (MDL), and maximum 
acceptability (MCL) are 0.20, 0.30, and 0.51, respectively, 
for maximum horizontal acceleration and 0.12. , 0.20 and 
0.34 were estimated for the maximum vertical acceleration 
and the maximum earthquake in the region with a 
magnitude of 7 (Abdan Faraz Consulting Engineers). In 
the study site, the earthquake coefficient has been 
determined and selected for stability analysis of 0.17 
(equivalent to one-third of the maximum. tolerated 
earthquake, based on Pyke's recommendation). Limiting 
the maximum acceleration of the input stimulus is 0.17 g 
due to the dynamic analysis of moderate and weak 
earthquakes by accepting low error and assuming linear 
soil behavior. It can also be related to the avoidance of 
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unrealistic tensile stresses in the lower shell elements, 
which occurred after adding dynamic to quasi-static 
stresses in earthquakes larger than 0.17 g due to the linear 
behavior assumption. Assuming the true nonlinear 
behavior of the soil and the inability to withstand the 
tensile strength of the aggregates, this problem can be 


solved. Therefore, to perform dynamic analysis and 
generate the input stimulus, the Tabas_ earthquake 
accelerometers with a maximum acceleration of 0.83 g and 
33 seconds have been used with the idea (Table 1 and 


Figure 3). 


Table 1. Earthquake characteristics used in dynamic analysis of Azadi Dam 
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Figure 3. Accelerometer used in dynamic analysis of Azadi Dam (Tabas earthquake) 


3. RESULTS AND DISCUSSION 


3.1, ANALYSIS OF STRESS AND STRAIN OF AZADI DAM IN STEADY-STATE 


The finite element analysis of Azadi Dam's body was 
performed after completing the construction phase and for 
steady-state conditions using Abaqus software. The model 
used in this analysis 1s Mohr-Columb elastoplastic, and the 


analysis is performed according to the parametric 
conditions of total stress. Also, the mechanical properties 
of Azadi Dam materials used in the modeling are shown in 
Table 2. 


Table 2. Geo mechanical parameters of stress-strain analysis 


ie — 

SS Fs, £f 

Zé AmMwm Bes 

sé PS~- S& 

= ay aS a. 
Core 0.72 ().42 
Filter 0.54 0.36 
Drain 0.52 0.34 
Shell 0.47 0.32 


Foundation 0.33 0.25 


Stress and strain modeling has been done by both quasi- 
static and dynamic methods in the full reservoir (steady- 
state). The quantities used in the analyzes are shown in 
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Table 3. The stress quantity is denoted by the symbol o, 
and the strain quantity is denoted by the symbol y. 
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Table 3. Symbol of stress and strain analysis committees of Azadi Dam in Abacus software 


Direction 


KK KKK x 


Sheet 


XZ 
XZ 
YZ 
XZ 
XZ 
YZ 


Quantity 


Oxx 
Oxy 
Oyy 
Yxx 
Yxy 
Yyy 


3.2. STRESS AND STRAIN ANALYSIS OF AZADI DAM BY THE QUASI-STATIC METHOD 


The simplest analysis of the behavior of a structure in 
an earthquake is the quasi-static method. This method 
is more common and older than other methods of 
seismic analysis. The use of this method dates back to 
1950. In fact, in this method, the effect of the 
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earthquake is considered statically in the analysis by 
applying forces that are obtained by multiplying the 
earthquake coefficients in the weight of the slippery 
mass along the horizon and vertical. Figure 4 shows the 
stress contour by the quasi-static method. 
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Figure 4. Distribution of Azadi dam stress in the quasi-static state in the direction of 6xx, Oxy and Gyy 


The results show that the highest stress in the oxx 
direction occurred due to the high adhesion of the clay 
and 483 kPa in the lower part of the core in the opposite 
direction of the X-axis. Also, the stress on both sides of 
the core (upstream and downstream shell) shows the 
value of 396 kPa, which is reduced by 18%, and as it 
moves away from the core, the amount of stress 
decreases to 135 kPa. The reduction of stress in this 
direction is due to heterogeneous materials. Most of the 
stress in the oxy state occurred in the upstream shell, 
which decreased by 32% with a certain geometric shape 
towards the upper levels. In this case, a lot of stress 





concentration has occurred and requires special 
measures; this amount of stress concentration in the 
upper shell is 122 kPa and in the lower shell is 154 kPa. 
The highest stress in the oyy state of 964 kPa occurred 
in both the upstream and downstream shells near the 
core. In the core, the stress of 386 kPa is created, which 
shows a 61% reduction compared to the shell. The 
reduction of stress in the shell and core upwards is not 
the same and the speed of this reduction in the core is 
higher than the shell. The strain contour of the Azadi 
dam in the quasi-static method is shown in Figure 5. 
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Figure 5. Distribution of Azadi dam strain in the quasi-static state in the direction of yxx, Yxy and yyy 


According to Figure 5, the greatest strain occurred in 
the yxx direction inside the core. The maximum strain in 
the core is due to fine-grained materials with low 
permeability, the value of which is equal to 0.0034, 
which increases towards the shell and reaches the value 
of 0.0002 on the embankment slope. The strain 
irregularity is greater in the y,y direction, due to the 
multiple behaviors of the heterogeneous materials of the 


dam, in this case, the maximum amount of strain in the 
core is 0.009 and the lowest amount of strain in the shell 
is 0.0015. The maximum strain in the yyy direction at the 
bottom of the core is 0.00132. The decrease in strain 
value of the yxx state is greater than the two states yxx 
and yxx. The reason for this behavior is the resistance of 
different parts of the heterogeneous dam. 


3.3. DYNAMIC STRESS AND STRAIN ANALYSIS OF AZADI DAM 


Dynamic analysis is the performance of numerical 
analysis on the body model and, if necessary, the 
foundation of the dam, which examines the behavior of 
the dam during the application of periodic seismic loads 
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and possible events after the earthquake by considering 
its behavior-strain of materials. The stress contour of 
the Azadi Dam is shown dynamically in Figure 6. 
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Figure 6. Distribution of Azadi dam stress in the dynamic state in the direction of xx, Oxy and Gyy 


The maximum amount of stress in the oxx direction after 
the seismic force in the core floor continues up to 30% of 
the shell width; this value is 502 kPa in the core and 640 
kPa in the shell, due to the coarse-grained materials of the 
dam shell and high hardness. This is the part. The amount 
of stress in the upstream and downstream shells increases 
towards the upper levels, decreasing at the point of 
collision with the core relative to the shell. The stress in 
the oxy direction occurs in the lower shell and 20% at the 
lower level of the core, which is between 171-145 kPa. 
Most stress changes in this condition occur in the lower 
Shell due to the phreatic line in this area. The core, in this 


8/ 


case, has almost regular stress of 42 kPa. The highest stress 
in the oyy direction occurs after the earthquake load 1s 
applied to the upstream shell, 15% away from the core, 
which is 1330 kPa. At the upper levels of the core at an 
angle of about 45°, the amount of stress is reduced so that 
in the middle of the shell will reach 650 kPa. Also, the 
stress in the lower shell and at the lower levels is the 
highest and decreases upwards. Its decreasing speed is 
higher than the upstream shell; the highest amount of stress 
in this part is 1059 kPa and in the middle of the shell is 
equal to 250 kPa. The strain contours of Azadi Dam are 
shown dynamically in Figure 7. 
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Figure 7. Distribution of Azadi dam strain in the quasi-static state in the direction of yxx, Yxy and Yyy 


In the case of strain in the direction of yxx, due to the 
fineness of the materials in this area and their high density, 
the highest strain occurred in the core after the earthquake, 
with the highest strain in this area being 0.0037 and the 
lowest strain being 0.0007. The strain, in this case, is 
almost the same in all upstream shells and is equal to 
0.0015, and the reason is the greater resistance of shell 
materials in this case. In the strain mode in the yxy 
direction, the maximum strain occurs as 0.0124 next to the 
core, then decreases to the inside of the core and reaches 
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the value of 0.003. In the upstream shell, the strain is 
approximately equal to 0.0011. The maximum strain 
occurs in the yyy direction at a distance of 0.125 of the 
width of the upstream shell next to the core and its value is 
approximately 0.016 and decreases with an angle of 
approximately 30 degrees to the higher level and reaches 
0.007 at the middle. The maximum strain in the core is 
0.014, which has decreased upwards. Three levels 1255m, 
1275m, and 1300m are used for shell stress diagrams 


(Figure 8). 
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Figure 8. Stress curve at different levels of Azadi Dam shell a) EL. 1300 m, b) EL. 1275 m and c) EL. 1255 m 


88 


J. Civil Eng. Mater.App. 2021 (June); 5(2): 81-92 


The interpretation of Figure 7 is shown in Table 4. 


Table 4. Stress values at different levels of Azadi Dam shell 
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610-640 2.14 413 4.6 910 6.6 1275 
145-190 1.8 50 9.8 410 6.8 1255 
According to Table 4, the highest stress occurs on the should be noted that in the first few seconds of an 
floor, which is 29% higher than the middle level and earthquake, the dam with a certain frequency (which is 
almost 2 times (1.85) compared to the upper level of the related to acceleration) shows resistance. Stress analysis in 
dam. The minimum stress value will be 72% higher than all three directions Gyy, Oxy, and oxx in the core is 
the middle and 4 times higher than the upper level. It investigated (Figures 10, 9, and 11). 
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Figure 9. Stress curve in 0x direction at different levels of Azadi Dam core a) EL. 1300 m, b) EL. 1275 m and c) EL. 
1255 m 
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Figure 10. Stress curve in oyy direction at different levels of Azadi Dam core a) EL. 1300 m, b) EL. 1275 m and c) EL. 
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Figure 11. Stress curve in oxy direction at different levels of Azadi Dam core a) EL. 1300 m, b) EL. 1275 m and c) EL. 
1255 m 


The interpretation of Figures 9,10 and 11 is shown in Table 5. 


Table 5. Stress values at different levels of Azadi dam core in the direction of oyy and Oxy, Oxx 
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As we expected, the maximum stress of the Azadi Dam 
core during the earthquake in all cases fell to low levels. 
The highest core stress is related to the oyy direction, which 
is 46% higher in the bottom, 49% in the middle, and 49% 
higher than the other directions. In all cases, the stresses 
increase and decrease with a limited frequency for 
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approximately 2.5 seconds, after which the stress values 
will reach their minimum and maximum with a relatively 
high slope. Strain diagrams are shown in Figures 12 to 14 
in all three directions yyy, Yxy and yxx at the core of the Azadi 
Dam. 


Figure 12. Strain curve in yxx direction at different levels of Azadi dam core a) EL. 1300 m, b) EL. 1275 m and c) EL. 
1255 m 
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Figure 13. Strain curve in yyy direction at different levels of Azadi dam core a) EL. 1300 m, b) EL. 1275 m and c) EL. 
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Figure 14. Strain curve in yx y direction at different levels of Azadi dam core a) EL. 1300 m, b) EL. 1275 m and c) 
EL. 1255 m 


The interpretation of Figures 12,13 and 14 is shown in Table 6. 


Table 6. Strain values at different levels of Azadi dam core in the direction of yyy and yxy, Yxx 
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The most strain occurs in the yxy direction and at high altitudes, almost twice as much as in other directions. 


4. CONCLUSION 

To control and prevent the increase of stresses during an 
earthquake, coarser-grained materials should be used in 
the body of the earth dam. The stress in the dynamic state 
is 49% higher in the ox direction, in the oxy 30% direction, 
and in the oyy 28% direction than in the quasi-static state. 
At the 1255m level, the maximum shell stress is 29%, at 
the 1275m level 68%, and at the 1300 m level 72% more 
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